Single-cell mtDNA heteroplasmy in colorectal cancer.

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

Optimal Management of Genetic Diversity in Subdivided Populations.

One of the main objectives of conservation programs is the maintenance of genetic diversity because this provides the adaptive potential of populations to face new environmental challenges. Genetic diversity is generally assessed by means of neutral molecular markers, and it is usually quantified by the expected heterozygosity under Hardy-Weinberg equilibrium and the number of alleles per locus or allelic diversity. These two measures of genetic diversity are complementary because whereas the former is directly related to genetic variance for quantitative traits and, therefore, to the short-term response to selection and adaptation, the latter is more sensitive to population bottlenecks and relates more strongly to the long-term capacity of adaptation of populations. In the context of structured populations undergoing conservation programs, it is important to decide the optimum management strategy to preserve as much genetic diversity as possible while avoiding inbreeding. Here we examine, through computer simulations, the consequences of choosing a conservation strategy based on maximizing either heterozygosity or allelic diversity of single-nucleotide polymorphism haplotypes in a subdivided population. Our results suggest that maximization of allelic diversity can be more efficient in maintaining the genetic diversity of subdivided populations than maximization of expected heterozygosity because the former maintains a larger number of alleles while making a better control of inbreeding. Thus, maximization of allelic diversity should be a recommended strategy in conservation programs for structured populations.

Gene Expression Analyses in Non Muscle Invasive Bladder Cancer Reveals a Role for Alternative Splicing and Tp53 Status.

Non-muscle invasive bladder cancer (NMIBC) represents a crucial problem for the national health care systems due to its high rates of recurrence and the consequent need of frequent follow-ups. Here, gene expression analyses in patients diagnosed as NMIBC were performed to determine those molecular pathways involved in tumor initiation, finding that both MYC and E2F are up regulated and helps to tumor initiation and progression. Our results also support an important involvement of alternative splicing events, modifying key pathways to favour bladder tumor evolution. Finally, since MDM2 showed differential exon usage, mutations in TP53 and its protein expression have been also studied in the same patients. Our data support that recurrence is epigenetically mediated and favoured by an increase protein expression of TP53, which appears more frequently mutated in advanced stages and grades, being associated to a worse prognosis. Therefore, TP53 mutational status could be used as a potential biomarker in the first stages of NMIBC to predict recurrence and prognosis.

metapop2: Re-implementation of software for the analysis and management of subdivided populations using gene and allelic diversity.

Management programmes often have to make decisions based on the analysis of the genetic properties and diversity of populations. Expected heterozygosity (or gene diversity) and population structure parameters are often used to make recommendations for conservation, such as avoidance of inbreeding or migration across subpopulations. Allelic diversity, however, can also provide complementary and useful information for conservation programmes, as it is highly sensitive to population bottlenecks, and is more related to long-term selection response than heterozygosity. Here we present a completely revised and updated re-implementation of the software metapop for the analysis of diversity in subdivided populations, as well as a tool for the management and dynamic estimation of optimal contributions in conservation programmes. This new update includes computation of allelic diversity for population analysis and management, as well as a simulation mode to forecast the consequences of taking different management strategies over time. Furthermore, the new implementation in C++ includes code optimization and improved memory usage, allowing for fast analysis of large data sets including single nucleotide polymorphism markers, as well as enhanced cross-software and cross-platform compatibility.

RNA-seq coupled to proteomic analysis reveals high sperm proteome variation between two closely related marine mussel species.

Speciation mechanisms in marine organisms have attracted great interest because of the apparent lack of substantial barriers to genetic exchange in marine ecosystems. Marine mussels of the Mytilus edulis species complex provide a good model to study mechanisms underlying species formation. They hybridise extensively at many localities and both pre- and postzygotic isolating mechanisms may be operating. Mussels have external fertilisation and sperm cells should show specific adaptations for survival and successful fertilisation. Sperm thus represent key targets in investigations of the molecular mechanisms underlying reproductive isolation. We undertook a deep transcriptome sequencing (RNA-seq) of mature male gonads and a 2DE/MS-based proteome analysis of sperm from Mytilus edulis and M. galloprovincialis raised in a common environment. We provide evidence of extensive expression differences between the two mussel species, and general agreement between the transcriptomic and proteomic results in the direction of expression differences between species. Differential expression is marked for mitochondrial genes and for those involved in spermatogenesis, sperm motility, sperm-egg interactions, the acrosome reaction, sperm capacitation, ATP reserves and ROS production. Proteins and their corresponding genes might thus be good targets in further genomic analysis of reproductive barriers between these closely related species. SIGNIFICANCE: Model systems for the study of fertilization include marine invertebrates with external fertilisation, such as abalones, sea urchins and mussels, because of the ease with which large quantities of gametes released into seawater can be collected after induced spawning. Unlike abalones and sea urchins, hybridisation has been reported between mussels of different Mytilus spp., which thus makes them very appealing for the study of reproductive isolation at both pre- and postzygotic levels. There is a lack of empirical proteomic studies on sperm samples comparing different Mytilus species, which could help to advance this study. A comparative analysis of sperm proteomes across different taxa may provide important insights into the fundamental molecular processes and mechanisms involved in reproductive isolation. It might also contribute to a better understanding of sperm function and of the adaptive evolution of sperm proteins in different taxa. There is now growing evidence from genomics studies that multiple protein complexes and many individual proteins might have important functions in sperm biology and the fertilisation process. From an applied perspective, the identification of sperm-specific proteins could also contribute to the improved understanding of fertility problems and as targets for fertility control.

Population genomics of parallel evolution in gene expression and gene sequence during ecological adaptation.

Natural selection often produces parallel phenotypic changes in response to a similar adaptive challenge. However, the extent to which parallel gene expression differences and genomic divergence underlie parallel phenotypic traits and whether they are decoupled or not remains largely unexplored. We performed a population genomic study of parallel ecological adaptation among replicate ecotype pairs of the rough periwinkle (Littorina saxatilis) at a regional geographical scale (NW Spain). We show that genomic changes underlying parallel phenotypic divergence followed a complex pattern of both repeatable differences and of differences unique to specific ecotype pairs, in which parallel changes in expression or sequence are restricted to a limited set of genes. Yet, the majority of divergent genes were divergent either for gene expression or coding sequence, but not for both simultaneously. Overall, our findings suggest that divergent selection significantly contributed to the process of parallel molecular differentiation among ecotype pairs, and that changes in expression and gene sequence underlying phenotypic divergence could, at least to a certain extent, be considered decoupled processes.

RNA-seq data from mature male gonads of marine mussels Mytilus edulis and M. galloprovincialis.

The mussels Mytilus edulis and Mytilus galloprovincialis are marine organisms with external fertilization able to hybridize where their distributions overlap allowing the study of reproductive isolation mechanisms in nature. We provide raw data of a transcriptomic analysis of mature male gonads from these two Mytilus spp. using NGS (Illumina) technology and a preliminary list of transcript that were functionally annotated showing species-specific differential expression. A shortlist including some of these genes and their corresponding proteins have been thoroughly analysed and discussed in Romero et al. (2018, Submitted for publication).

On the Consequences of Purging and Linkage on Fitness and Genetic Diversity.

Using computer simulation we explore the consequences of linkage on the inbreeding load of an equilibrium population, and on the efficiency of purging and the loss of genetic diversity after a reduction in population size. We find that linkage tends to cause increased inbreeding load due to the build up of coupling groups of (partially) recessive deleterious alleles. It also induces associative overdominance at neutral sites but rarely causes increased neutral genetic diversity in equilibrium populations. After a reduction in population size, linkage can cause some delay both for the expression of the inbreeding load and the corresponding purging. However, reasonable predictions can be obtained for the evolution of fitness under inbreeding and purging by using empirical estimates of the inbreeding depression rate. Purging selection against homozygotes for deleterious alleles affects the population's pedigree. Furthermore, it can slow the loss of genetic diversity compared to that expected from the variance of gametic contributions to the breeding group and even from pedigree inbreeding. Under some conditions, this can lead to a smaller loss of genetic diversity, even below that expected from population size in the absence of selection.

Allelic diversity for neutral markers retains a higher adaptive potential for quantitative traits than expected heterozygosity.

The adaptive potential of a population depends on the amount of additive genetic variance for quantitative traits of evolutionary importance. This variance is a direct function of the expected frequency of heterozygotes for the loci which affect the trait (QTL). It has been argued, but not demonstrated experimentally, that long-term response to selection is more dependent on QTL allelic diversity than on QTL heterozygosity. Conservation programmes, aimed at preserving this variation, usually rely on neutral markers rather than on quantitative traits for making decisions on management. Here, we address, both through simulation analyses and experimental studies with Drosophila melanogaster, the question of whether allelic diversity for neutral markers is a better indicator of a high adaptive potential than expected heterozygosity. In both experimental and simulation studies, we established synthetic populations for which either heterozygosity or allelic diversity was maximized using information from QTL (simulations) or unlinked neutral markers (simulations and experiment). The synthetic populations were selected for the quantitative trait to evaluate the evolutionary potential provided by the two optimization methods. Our results show that maximizing the number of alleles of a low number of markers implies higher responses to selection than maximizing their heterozygosity.

The action of stabilizing selection, mutation, and drift on epistatic quantitative traits.

For a quantitative trait under stabilizing selection, the effect of epistasis on its genetic architecture and on the changes of genetic variance caused by bottlenecking were investigated using theory and simulation. Assuming empirical estimates of the rate and effects of mutations and the intensity of selection, we assessed the impact of two-locus epistasis (synergistic/antagonistic) among linked or unlinked loci on the distribution of effects and frequencies of segregating loci in populations at the mutation-selection-drift balance. Strong pervasive epistasis did not modify substantially the genetic properties of the trait and, therefore, the most likely explanation for the low amount of variation usually accounted by the loci detected in genome-wide association analyses is that many causal loci will pass undetected. We investigated the impact of epistasis on the changes in genetic variance components when large populations were subjected to successive bottlenecks of different sizes, considering the action of genetic drift, operating singly (D), or jointly with mutation (MD) and selection (MSD). An initial increase of the different components of the genetic variance, as well as a dramatic acceleration of the between-line divergence, were always associated with synergistic epistasis but were strongly constrained by selection.

Analysis and management of gene and allelic diversity in subdivided populations using the software program METAPOP.

METAPOP (http://webs.uvigo.es/anpefi/metapop/) is a desktop application that provides an analysis of gene and allelic diversity in subdivided populations from molecular genotype or coancestry data as well as a tool for the management of genetic diversity in conservation programs. A partition of gene and allelic diversity is made within and between subpopulations, in order to assess the contribution of each subpopulation to global diversity for descriptive population genetics or conservation purposes. In the context of management of subdivided populations in in situ conservation programs, the software also determines the optimal contributions (i.e., number of offspring) of each individual, the number of migrants, and the particular subpopulations involved in the exchange of individuals in order to maintain the largest level of gene diversity in the whole population with a desired control in the rate of inbreeding. The partition of gene and allelic diversity within and between subpopulations is illustrated with microsatellite and SNP data from human populations.

Methylation changes associated with early maturation stages in the Atlantic salmon.

BACKGROUND: Early maturation in the Atlantic salmon is an interesting subject for numerous research lines. Prior to sea migration, parr can reach sexual maturation and successfully fertilize adult female eggs during the reproductive season. These individuals are known as precocious parr, mature parr or "sneakers". Reasons for early maturation are unknown and this transitory stage is usually considered to be a threshold trait. Here, we compare methylation patterns between mature and immature salmon parr from two different rivers in order to infer if such methylation differences may be related to their maturation condition. First we analyzed genetic differences between rivers by means of AFLPs. Then, we compared the DNA methylation differences between mature and immature parrs, using a Methylation-Sensitive Amplified Polymorphism (MSAP), which is a modification of the AFLPs method by making use of the differential sensitivity of a pair of restriction enzymes isoschizomeres to cytosine methylation. The tissues essayed included brain, liver and gonads. RESULTS: AFLPs statistical analysis showed that there was no significant differentiation between rivers or a significant differentiation between maturation states in each river. MSAP statistical analysis showed that among the three tissues sampled, the gonads had the highest number of significant single-locus variation among populations with 74 loci followed by brain with 70 and finally liver with only 12. Principal components analysis (PCA) of the MSAP profiles revealed different profiles among different tissues (liver, brain and testis) clearly separating maturation states in the testis tissue when compared to the liver. CONCLUSIONS: Our results reveal that genetically-similar mature and immature salmon parr present high levels of DNA methylation variation in two of the three analyzed tissues. We hypothesize that early maturation may be mostly mediated by epigenetic processes rather than by genetic differences between parrs. To our knowledge this is the first study that attempt to link phenotypic plasticity in salmonids and epigenetic changes.

Early detection of population declines: high power of genetic monitoring using effective population size estimators.

Early detection of population declines is essential to prevent extinctions and to ensure sustainable harvest. We evaluated the performance of two N e estimators to detect population declines: the two-sample temporal method and a one-sample method based on linkage disequilibrium (LD). We used simulated data representing a wide range of population sizes, sample sizes and number of loci. Both methods usually detect a population decline only one generation after it occurs if N e drops to less than approximately 100, and 40 microsatellite loci and 50 individuals are sampled. However, the LD method often out performed the temporal method by allowing earlier detection of less severe population declines (N e approximately 200). Power for early detection increased more rapidly with the number of individuals sampled than with the number of loci genotyped, primarily for the LD method. The number of samples available is therefore an important criterion when choosing between the LD and temporal methods. We provide guidelines regarding design of studies targeted at monitoring for population declines. We also report that 40 single nucleotide polymorphism (SNP) markers give slightly lower precision than 10 microsatellite markers. Our results suggest that conservation management and monitoring strategies can reliably use genetic based methods for early detection of population declines.

The action of purifying selection, mutation and drift on fitness epistatic systems.

For different fitness mutational models, with epistasis introduced, we simulated the consequences of drift (D scenario) or mutation, selection, and drift (MSD scenario) in populations at the MSD balance subsequently subjected to bottlenecks of size N = 2, 10, 50 during 100 generations. No "conversion" of nonadditive into additive variance was observed, all components of the fitness genetic variance initially increasing with the inbreeding coefficient F and subsequently decreasing to zero (D) or to an equilibrium value (MSD). In the D scenario, epistasis had no appreciable effect on inbreeding depression and that on the temporal change of variance components was relevant only for high rates of strong epistatic mutation. In parallel, between-line differentiation in mean fitness accelerated with F and that in additive variance reached a maximum at F approximately 0.6-0.7, both processes being intensified by strong epistasis. In the MSD scenario, however, the increase in additive variance was smaller, as it was used by selection to purge inbreeding depression (N > or = 10), and selection prevented between-line differentiation. Epistasis, either synergistic or antagonistic (this leading to multiple adaptive peaks), had no appreciable effect on MSD results nor, therefore, on the evolutionary rate of fitness change.

Preserving population allele frequencies in ex situ conservation programs.

Optimization of contributions of parents to progeny by minimizing the average coancestry of the progeny is an effective strategy for maintaining genetic diversity in ex situ conservation programs, but its application on the basis of molecular markers has the negative collateral effect of homogenizing the allelic frequencies at each locus. Because one of the objectives of a conservation program is to preserve the genetic composition of the original endangered population, we devised a method in which markers are used to maintain the allele frequency distribution at each locus as closely as possible to that of the native population. Contributions of parents were obtained so as to minimize changes in allele frequency for a set of molecular markers in a population of reduced size. We used computer simulations, under a range of scenarios, to assess the effectiveness of the method in comparison with methods in which contributions of minimum coancestry are sought, either making use of molecular markers or genealogical information. Our simulations indicated that the proposed method effectively maintained the original distribution of allele frequencies, particularly under strong linkage, and maintained acceptable levels of genetic diversity in the population. Nevertheless, contributions of minimum coancestry determined from pedigree information but ignoring the genealogy previous to the conservation program, was the most effective method for maintaining allelic frequencies in realistic situations.